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ABSTRACT 

Various balloon and satellite observatories have revealed what appears to be 
an extended source of 0.511 MeV annihilation radiation with flux of ~ 10 -3 
photons cm -2 s -1 centered on the Galactic Center. Positrons from radioactive 
products of stellar explosions can account for a significant fraction of the emission. 

We discuss an additional source for this emission: namely e + e~ pairs produced 
when X-rays generated from the ~ 2.6 x 10 6 solar mass Galactic Center Black 
Hole interact with ~ 10 MeV temperature blackbody emission from 10 17 g black 
holes within 10 14-15 cm of the center. The number of such Small Mass Black 
Holes (SMMBHs) can account for the production of the IQ 42 e + s -1 that produces 
the observed annihilation in the inner Galaxy when transport effects are taken 
into account. We consider the possibility for confirming the presence of these 
SMMBHs in the Galactic Center region with future generations of gamma-ray 
instruments if a blackbody like emission of ~ 10 MeV temperature would be 
detected by them. Small Mass Black Hole can be a potential candidate for dark 
(invisible) matter halo. 
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1. Introduction 

The Galactic Center (GC) of our Galaxy was detected in 1974 as a strong radio source 
called Sgr A* and is the site where a supermassive black hole with 2.6 xlO 6 M© is present 
(e.g. Melia and Falcke 2001). Morever the total mass of stars enclosing the inner 20 pc of the 
Galaxy contains ~ 10 8 M©, or ~ 0.1% of the Galactic stellar mass (Baganoff et al. 2003). 
Therefore, the X-ray observations of the Galactic Center can potentially reveal accreting 
black holes and neutron stars. During five years of Chandra observations, Baganoff et al. 
identified seven X-ray transients located within 20 pc of Sgr A*. Remarkably four of these 
transients lie within only 1 pc of Sgr A* (Muno et al. 2005). This implies that, compared to 
the numbers of similar systems located between 1 and 20 pc, transients are over-abundant 
by a factor of ~ 20 per unit stellar mass within 1 pc of Sgr A*. 

Studies of the GC are also a subject of intense activity due to the possible presence 
of a dark matter cusp. Numerous space missions have observed the center of our Galaxy 
in radio, infrared, X and 7 rays. Thanks to these measurements we have a more precise 
picture of this region. There is a strong activity both in X and 7 rays near the center of our 
Galaxy. A recent HESS (High Energy Stereoscopic System) experiment reported a detection 
of very high energy gamma rays from the Galactic center within 1' of Sgr A*. The flux above 
165 GeV is order of 1.8 x 10~ 7 mr 2 s"" 1 and consistent with a luminosity of 10 35 ergs s -1 
(Aharonian, et al. 2004). 

The most intense galactic gamma-ray line is 511 keV annihilation line. Balloon-flight 
observations first detected the annihilation line from the Galactic Center region (Johnson, 
Harnden & Haymes 1972; see also Johnson k, Haymes, 1973). Since that time the inner 
Galaxy has been the object of intense observations from space missions and balloon ex- 
periments (Leventhal, MacCallum & Stang 1978; Purcell et al. 1993; Ramaty, Skibo & 
Lingenfelter 1994; Purcell et al. 1997; and Cheng et al. 1997; Share et al. 1990; Teegarden 
et al. 1996; Gehrels et al. 1991; Leventhal et al. 1993; Niel et al. 1990; Chapuis et al. 1991). 

Recently the new measurements of the 511 keV line emission from the Galactic Center 
region have been performed with spectrometer SPI on the space observatory INTEGRAL 
(Jean et al. 2003). The INTEGRAL observation is in a good agreement with the previous 
measurements and is well fitted by a gaussian with the full width at half maximum of ~ 
10 degrees, with a 2 a confidence interval of 6-18 degrees. The INTEGRAL measurement 
providing the line width of ~ 2.4 keV exlcudes velocity v of the interstellar medium (ISM) 
greater than ~ 800 km s _1 , i.e. /3 = v/c < 2 x 10 -3 (Churazov et al. 2005). The source line 
flux is : 9.9 x 10" 4 ph cm -2 s _1 , for which the corresponding luminosity of the line emission 
is about 6 x 10 36 ergs s _1 for the distance of about 8 kpc. Then the positron sources are 
produced at a rate given by; ~ 10 42 d| s -1 , where d 8 is the distance expressed in unit of 8 
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kpc. In our calculations (see §4) we use this value of the pair flux 10 42 s 1 inferred from the 
INTEGRAL observations. 

Several possible sources of the observed line related to the positron sources have been 
proposed. Most of these propose that positrons come from radioactive elements ejected by 
stellar explosions and massive stars (Matz et al. 1988; Kurfess et al. 1992, Diehl et al. 1995; 
Iyudin et al. 1994; Dermer & Skibo 1997). Boehm, Fayet &; Silk (2004) offer an alternative 
explanation that the observed annihilation line could be related to the annihilation of dark 
light matter particles into e + e _ . Guessoum, Jean & Gillard (2005) review the main processes 
that lead to the positron production. They are: (i) the fJ + decay of radioactive nuclei which 
is widespread in the explosive and/or hydrostatic nucleosynthesis enviroments of novae, 
supernovae, Wolf-Rayet, and Asymptotic Giant Branch stars, (ii) 7r+ decay into p, + and 
positrons in collisions of highly energetic ( E > 200 MeV) cosmic rays with interstellar 
material, (iii) pair production by interaction of an electron with a strong magnetic field and 
(iv) photon-photon pair production of high energy photons near luminous compact objects, 
black holes, active nuclei (AGN) etc. The latter origin is consistent with production from 
interaction of X-ray and 7— emission in the very compact shell near Galactic Center as we 
discuss in this paper. 

X-ray observations of the Galactic Center have been made with several instruments 
(Watson et al. 1981; Skinner et al. 1987; Kawai et al. 1988; Takeshima et al. 1991; Sunyaev 
et al. 1993; Markevitch et al. 1993; Pavlinsky et al. 1993; Preddehl, & Trumper 1994). High- 
sensitivity imaging observations with the Einstein Observatory (Watson et al.) and ROSAT 
(Preddehl & Trumper) resolved several discrete sources within 10' of Galactic Center. One 
of the ROSAT sources, RXJ 1745.6-2900, is within 10" coincident of the nonthermal radio 
source Sgr A *, which is believed to be the dynamical center of our Galaxy. 

The Japanese satellite GINGA detected K a lines from highly ionized iron, suggesting 
a large amount of hot plasma strongly concentrated towards the Galactic Center (Koyama 
et al. 1989; Yamauchi et al. 1990). The ASCA X-ray image (Tanaka et al. 1995) has 
revealed extended diffuse emission over the Galactic Center region. The observed spectrum 
shows the K Q lines from highly-ionized ions of silicon, sulphur, argon, calcium, in addition 
to high-ionization iron lines. This gives firm evidence for the presence of high-temperature 
(kT ss 10 keV) plasma. A large energy input of ~ 10 41 “ 42 ergs s _1 is required to produce 
the observed plasma. Koyama et al. (1996) argue that the source of this energy resides at 
the Galactic Center. The luminosity required to account for the 6.6 keV line from Sgr B2 
(the prominent molecular cloud) is 2 x 10 39 ergs s _1 (in 2-10 keV band). While the Sgr A 
region has an X-ray luminosity of ~ 10 36 ergs s” 1 at present, it could have been much more 
luminous ~ 300 yr ago. The distinctly bright region inside the Sgr A shell is probably due 
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to this recent activity and the X-ray luminosity may have been as high as several xlO 39 
ergs s _1 at early times. Koyama et al. conclude that the ASCA results suggest the presence 
of an active galactic nucleus and a large mass concentration at the Galactic Center with 
transient activities due to sporadic mass accretion rate. Muno et al. (2004a) also presented 
the study of the spectrum of diffuse emission in IT x 17' field around Sgr A* during 6 x10 s 
s of Chandra observations. They confirmed and extended the results of the previous ASCA 
study. In particular, they showed that the continuum and line spectra are consistent with an 
origin in a two-temperature plasma: As a cooler, kT & 0.8 keV plasma with varying surface 
brightness, and a hotter, kT « 8 keV component that is more spatially uniform. 

Recently Revnivtsev et al. (2004) have published data on fluorescence in giant molecular 
cloud Sgr B2’s (ICR J17475-2822). They have constructed the broad band (3-200 keV) 
spectrum from the ASCA-GRANAT-INTEGRAL data. It appears that the X-ray irradiating 
spectrum of the cloud is, in fact, not free-free, (as early proposed above) but it is rather quite 
hard, extending to over 100 keV with a power law with photon index 1.8. Revnivtsev et 
al. concluded that 300-400 years ago Sgr A* was a low luminosity (of ~ 2 x 10 39 ergs s _1 
at 2-200 keV) AGN with characteristic hard X-ray Comptonization spectrum with photon 
index of about 1.8 and plasma temperature of about 50 keV. 

Gamma-ray emission between 1 MeV and 30 MeV and higher has been measured by 
COMPTEL and EGRET respectively within a few degrees of the GC cusp. Using the 
COMPTEL (Strong et al. 1998, 1999, 2004, 2006) and EGRET data (Hunter et al. 1997) 
one can calculate that the luminosity in the 1-100 MeV energy band is about 4 x 10 38 ergs 
s -1 assuming that the distance to the Galactic Center is about 8 kpc. In the energy band 
between 1 and 30 MeV the COMPTEL spectrum presented in units of E 2 photon intensity 
cm -2 sr -1 s -1 MeV is quite flat (taking into account the error bars). The representative 
value of the flux in this energy band is about 2 x 10~ 2 cm -2 sr -1 s -1 MeV. 

In this Paper, we study pair creation due to X-ray and 7*- photon interactions inside the 
compact Galactic Center shell (the SGR A shell) and we discuss the possibility of annihilation 
line formation when these pairs propagate through the Galaxy from the Galactic Center 
region. In §2 we present a pair creation model. In §3 and §4 we compare our theoretical 
predictions for the pair and annihilation photon production with observations. In §3 we 
determine the size of the pair production area using the observable fluxes in X-ray and 
7— ray energy ranges. In §4 we investigate the outward propagation of positrons from the 
Galactic Center region. In §5 we present arguments for the origin of the ~ 10 MeV emission 
near Galactic Center. We discuss our results in §6 and we make our final conclusions in §7. 



2. Pair creation model 


In Figure 1 we present our model for the illumination of the 1-100 MeV 7— ray point like 
sources by X-ray radiation of the central black hole. We assume that the photon density of 
X-ray photons emitted by the inner shell near the central BH within 10 3 Rs (Schwarzschild 
radii) is described by the following distribution over radius r 

Ux = 4tt r*EhTc ^ 

where Ex is X-ray luminosity of the central source per keV (X-ray energy spectrum), Ex is 
the energy of X-ray photon and c is the speed of light. Also we consider that the spherical 
shell between the inner radius R c and R (R c < r < R and R R c ) is uniformly filled by 
7— radiation sources where photon density production P 7 is 


E^4irR 3 /3 ^ 

where £ 7 is the 7-radiation luminosity of the spherical shell per keV (7-ray energy spectrum) 
and E 7 is the energy of the 7-ray photon. Then the pair production by the photons (for 
which the pair production condition ExE-y > ( m e c 2 ) 2 is satisfied) is 


P X J] + = 4t r [ R r 2 P-ydr [ T* )7 (r, tp, 6)dQ (3) 

J R c J Q 

where Tx tl {r, <p, 6) is the optical path for 7-ray photon propagating in the X-ray background 
and incoming at radius r at given direction O. The multiple integral in equation (3) can be 
calculated analytically. In fact, 


f T Xn (r, tp, 6)dQ = {2ttox 7 } x 
Jn 

f p 1 I" poo pr 1 poo ^ 

< dp, n(r')ds + / n(r')ds + / dp, / n(r')ds> (4) 

WO .J r(l“/i 2 ) 1 / 2 Jr(l— fi 2 ) 1 / 2 _ J — 1 Jr J 


where 


n(r')ds 


n^r'dr' 


( 5 ) 


r' 2 [r' 2 — r 2 ( 1 — p 2 )] 1 / 2 

and n 0 — Cx/4xExC] where ax,^{Ex, E^) is a pair creation cross section given by (Dirac 
1930, Heitler, 1954) 

Kill 


vx, 1 {E x ,E 1 ) = — T^(l - b 2 ) 


2 b{b 2 - 2) + (3 - 6 4 ) In 


1 + b 
1-b 


2 


( 6 ) 
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and b = [1 — (m e c 2 ) 2 / ExE^] 1 ? 2 , and where the classical electron radius TZa — e 2 /ra e c 2 
2.818 x 10" 13 cm. 

The analytical integration of f n{r')ds 

dr' 


J n{r')ds =5 no J — 


' \/r'^ — r 2 (l — /i 2 ) r-y/I — /r 2 
allows us to calculate the second integral of f n(r')ds over /n 

Lx 3 tt 2 


n 0 r^/l -yu 2 

arccos y , — — h C 


[ Tx tl (r, ip, 0)dQ = 27R7x, 7 - 


/n A-kExc 4r 

Finally, the expression for the Green’s function P*’J + can be written as 

pX , 7 97 t 2 ^Xj^LxLj 

e ~ e+ ~ 1 6cR E X E^ ' 


(7) 


( 8 ) 


(9) 


In order to obtain the total pair production rate P e - e + one should integrate P*’J + over X-ray 
and 7— energies, namely 

9tt 2 /, °° 


1 V 


e~e+ 


16cP 


[°°dEx [ 

J 0 «/{? 


(m e c 2 ) 2 /Kjt 


Ex Ex 7 * * * * 


( 10 ) 


3. Results of calculations 

To compute the positron flux one needs the X and 7 ray energy spectra, Lx and L 1 

respectively. We will use for the total X-ray and 7 luminosities of the inner GC region, 

Lx ~ 2 x 10 39 ergs s _1 (Revnivtsev et al. 2004) and L 7 ~ 4 x 10 38 ergs s -1 (e.g. Strong et 

al. 1999) respectively (see a review of the observations in the Introduction). 

We assume that the 7— radiation emanates from a optically thick medium (as an ac- 

cretion disk around the compact object) and consequently the emergent spectrum has a 

blackbody like shape that color temperature T 7 = 10 MeV (see Fig. 1). Then for a normal- 

ization constant C 1 we have an equation as follows 


L*y — f L-y(E-y)dEj — Cj f 
Jo Jo 


00 y 3 


exp (y) 


-dy 


( 11 ) 


where y = E^/kT. ' T For the X-ray photons, we assume a power law with an exponential 
cutoff 

poc poo 

Lx = Lx{Ex)dEx — Cx / x~ a exp(— x)dx (12) 

Jo Jo 
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where x — E x jkT x , kT x = 100 keV, the spectral index a = l' — 1 = 0.8 and C x is the 
normalization constant. From the previous section we saw that the process of pair creation is 
driven by the cross section given by equation (6). The energy dependence, represented by the 
parameter b can be expressed using the dimesionless energies x = E x / kT x and y — Ej/kTy, 
namely b=[l — (m e c 2 ) 2 / E X E^ 2 = [1 — 3.26/xy] 1 ^ 2 . Then the formula for the pair flux 


e~e+ — 


j* i r dE r ^(Ex,£ 7 )£ 7 (c,)^(Sx) Jp 

IteRJo X hv*e-M* ExE, 


can be modified using the variables x and y as follows 


P, j-e+ — 


f°° dx f°° <'B x ^,(x,v)C,(y)C x M j 
16 CRJ 0 A. 26/1 X V 


( 13 ) 


(14) 


We compute this double integral and obtain the following result: 

3.2 x 10 56 


e~e+ 


R 


cm 1 s 1 . 


(15) 


To reproduce the value of 10 42 s 1 for the positrons rate, one needs a value of the radius R, 
given by: 


R = 


3.2 x 10 56 * 
Pe~e+ 


3.2 X 10 56 
lO 42 * 


- 3.2 x 10 14 


cm. 


(16) 


It is worth noting that the intrinsic pair flux ~ 10 36 s _1 , formed as a result of 7 + 
7 interaction (to calculate this flux one should replace by £ 7 in formula 14) in this 
7— emission region of size 3.2 x 10 14 cm, is six orders of magnitude less than the pair flux 
~ 10 42 s _1 required to reproduce the observed annihilation line strength. 


4. Outward propagation of positrons from the Galactic Center region 

A positron, e + may collide with an electron e~, to produce two gamma ray photons 

according to the reaction e + + e~ — > 7 + 7. Before decaying, the positron in free space can 

interact with an electron and to form positronium, a bound state with lifetime of order of 

10 -10 s. The particles annihilate into two or three photons depending whether it is formed 

in the singlet or triplet state. In our calculations we assume that most of the positrons are 

thermalized in the surrounding thermal plasma and thus we can consider the annihilation of 

thermal positrons at a temperature of order 10 7 K , i.e. =< v > jc ^ 10 -3 , where < v > 

is the thermal velocity of electrons. In fact, the recent INTEGRAL observations of 511 keV 

annihilation line (Churazov et al. 2005) exclude velocities of the interstellar medium greater 
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than 10 3 c. For /3 «C 1, the cross-section a, for positron annihilation with a free electron at 
rest is given by Dirac (1930) 

cr « TKrlifi. (17) 

The problem of the positron thermalization and annihilation photon production in the in- 
terstellar medium has been recently studied in details by Guessoum, Jean, & Gillard (2005) 
(see for a review of this subject). They reexamine in detail the various processes for positron 
annihilation using most recent interaction cross section with atomic and molecular hydrogen, 
as well as helium. Putting all the new calculations together, they constructed annihilation 
spectra of ISM. 

The fundamental question is how far the positrons of energy about one MeV can travel 
from the Galactic central region to be annihilated. In order to answer to this question one 
has to compare the diffusion (travel) time with the energy loss (thermalization) time taking 
into account the main energy loss mechanisms. In general terms, the particle diffusion is 
related to escape of a magnetized turbulent plasma. 

Jean et al. (2005) analyze the positron diffusion propagation in the GC region in detail. 
The quasilinear diffusion coefficient D q i (Melrose 1980) can be expressed as: 

D„ = D b (^-) r' (18) 

\ Amax J 


with Dq = (1/3) Wl the Bohm diffusion coefficient, X max is the maximum scale of the 
turbulence, 5 = 5/3 for a Kolmogorov turbulent spectrum, rj — 5B 2 / < B 2 > is relative 
perturbation of the magnetic field pressure which is often approximated to 1, v is the positron 
velocity and 77 is the Larmor radius (gyroradius) 


r L = 


m e c 'yv 


(19) 


q is the positron (electron) charge, B is the magnetic field strength, 7 = [1 — (v/c) 2 ] -1 / 2 . 
The maximum scale X max was estimated to be ~ 100 pc from the measurements of ISM 
turbulence (Armstrong et al. 1995). 


Recently the (mean) magnetic field strength in the GC region was measured by La Rosa 
et al. (2005), hereafter LaR05, to be ~ 10 /xG. LaR05 have used the Very Large Array in all 
four configurations to image the GC region at 74 MHz. The resulting image of large scale 
(6° x 2°) of nonthermal synchrotron emission presented in LaR05 has a resolution of 125". 
At distance of 8 kpc the angular scales of 6° x 2° corresponds to region 840 pc x 280 pc 
and that of 125" corresponds to 0.09 pc. It has be noted that in the past several authors 
(see Morris & Serabyn et al. 1996 and references therein) estimated magnetic fields values 
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of ~ mG in the GC region. This is 2 orders of magnitude larger than the value (~ 10 fj,G) 
obtained by LaR05. 

Thus for the mean B of 10 /iG and for 1 MeV positrons D q i ~ 2.6 x 10 26 cm 2 s -1 . The 
low limit time scale for 1 MeV positrons to thermalize in the B ~ 10 /j,G and ~ 1 cm " 3 
density region due to ionization, bremsstrahlung, synchrotron and other processes is r> 10 s 
years. Then the distance travelled by diffusion in a time of ~ 10 7 years (slowing down 
time plus annihilation time in the standard grain model) is 

dqi — yj DtfTd - 3 kpc (20) 

which is much greater than typical half-size of the hot GC region. The volume filling factor 
of the hot GC region 0.72 is highest among all phases of the Galactic bulge (see more details 
of these estimates in Jean et al. 2005). Using the spectral data and positron propagation 
analysis Jean et al. (2005) come to conclusion that the sources of annihilation emission 
are diffusively distributed over the Galaxy. They also explain that the lack of annihilation 
emission from the GC hot gas is due to its low density, which allows positrons to escape this 
phase. 

It is worth noting Liu, Petrosian & Melia (2004) find that an outflow of low energy 
electrons of order 1 MeV are distributed over a spatially larger scale than that of the BH 
inner region (see more discussion of their work in section 6.2). Namely, the escape time of 
these electrons is more than a factor of 2 shorter than the time scales for acceleration and 
losses. 

On the hand one can argue that the positrons can be diffusely trapped, thermalized and 
annihilated in a region of size L of order a few parsec where the magnetic field can be much 
higher than ~ 10 jiG. Melia & Falcke (2001) suggest that out to ~ 2 — 3 pc (essentially 
inside the cavity surrounded by the circumnuclear disk) , the field is turbulent (matching the 
turbulent plasma generated by wind-wind interactions from the Wolf-Rayet and O/B stars 
in this region), and can be as high as ~ 0.2 — 1 mG. 

The problem with such fields is that the gyroradius of the electrons and positrons 
produced in gamma-gamma interactions, will be far smaller than the region over which 
the positrons are believed to be annihilating. Consequently, the pairs should be radiating 
profusely via synchrotron, bremsstrahlung, and inverse Compton. 

Another estimate of the magnetic field comes from equipartition with the photon density. 
The photon energy density in this region is about 1 eV cm -3 . To obtain this number we 
assume that the bolometric luminosity of the source in the Galactic Center is of order 10 38 
erg s -1 (see e.g. Narayan et al. 1998). Thus the derived magnetic strength i 3 bo i using 
equipartition between photons and magnetic field energies is about 10 /iG. For such a low 
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magnetic field we have already demonstrated that the positrons can escape from the GC 
region. 

Longair (1994) also derived the minimum magnetic energy requirement for a given lumi- 
nosity C u and photon volume V per unit time (in our case V = 4 ttL 2 c). He shows that B m i n oc 
(C. v /V) 2 ! 5 * 7 a value that is higher than that obtained using the relation £^/(47 tL 2 c) = B 2 ol /8n, 
But the ratio of B min j B M is not large. The indicies of the ratios C u /V and C v /{A / nL 2 c) are 
2/7 to 1/2 respectively. Furthermore Longair found that the magnetic field strength, B m i n , 
occurs when there is near equality of the energies in the particles and magnetic field. LaR05 
confirm this conclusion. They show that the mean magnetic field in the GC region 10 //G 
inferred from their observations is consistent with the particle energy density of 1,2 eV found 
in cosmic-ray data. 

The turbulent diffusion time scale r — L 2 lD ql is not only determined by the magnetic 
field. In fact for 5 — 5/3 the diffusion coefficient D ql is proportional to B~ 1 ^ (see Eq.18). 
Even if the magnetic field is about 1 mG the diffusion coefficient D q i decreases by factor 4.6 
only, namely D q i ~ 5.5 x 10 25 cm 2 s -1 , The related difussion time through this region r is 
about 5 x 10 4 years which is much less than the thermalization and annihilation time 10 7 
years. 


5. Emission from small mass BHs as a source of 7 — radiation 

One of the possible origins for MeV-emission is a disk emission from mass accretion by 

Small Mass Black Holes (SMMBHs) . The BH mass can be evaluated using the color disk 
temperature T coi . Using formula (5) in Shrader & Titarchuk (1999), hereafter ShT99, we 
obtain that a black hole mass in solar units (m — M^/Mq) is; 


rn = m 


(i/7)3 4 r 4 

(Tmax) 4 [(7 /6) 2 r j n ] 3 


( 21 ) 


77 i = M/M Edd is the dimensionless mass accretion rate in units of M Edd == L Ed d /c 2 , related 
to the Eddington luminosty L^dd, T max k, 1.2 T co i is the maximum temperature in the disk 
in keV, Th — T co i/T e ff is a disk color factor and r in is the inner disk radius calculated in 
Schwarzchild units. Analyzing quite a few BH sources, ShT99 and Borozdin et al. (1999), 
hereafter BRT99, established that in the soft state the inner disk radius is very close to the 
last stable Keplerian orbit, namely r in « 3. The disk color factor Th « 2.6 was calculated 
by BRT99 using the known contraints on BH mass and the distance to the source for GRO 
J1655-40. This value of the disk color factor Th ~ 2.6 obtained for GRO J1655-40 was 
recently confirmed for other Galatic black hole candidates sources. Shrader & Titarchuk 
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(2003) have made use of improved Galactic black hole binary dynamical mass determinations 
to derive, in effect, an empirical calibration of this factor. 


In the soft state the mass accretion rate in the disk m = M / M Edt j is on the order of 
Eddington, i.e m«l. With assumption regarding rj n , Th and rh we obtain 

(l/7)3 4 (T ft /2.6) 4 


m = 3.7 x 10 16 (rh/l) 


CW 1.2 • 10 4 keV) 4 [(7/6) 2 (rj„/3] 3 ’ 


( 22 ) 


The SMMBHs can be exposed through the accretion, if the mass accretion is at the 
Eddington rate and higher. The Eddington accretion rate by definition is 

Mem = L E dd/c 2 = 1.4 x 10 17 m g s -1 . (23) 

Thus for SMMBH of m = rn STnm bh = 3.7 x 10 -16 (the SMMBH mass in grams is ~ 7 x 10 17 
g) the Eddington rate is about 52 g s -1 . The gravitation attraction of the black hole in 
the Galactic center causes a substantial mass accretion in its proximity. This accretion flow 
passing through SMMBH which are probably distributed uniformly around the central black 
hole, provides enough material for mass supply of the SMMBH objects. 

Sgr A*, the massive BH in the Galactic Center was recently found to be surrounded 
by a cluster of young, massive stars. Over the past year, three of these stars named SO-2, 
SO-16 and SO-19 by Ghez et al (2004) and Schodel et al (2003) were determined to be at 
~ 10 3 Schwarzschild radii from the central black hole. It was also found that these stars can 
be classified as O stars in the main sequence (Ghez et al. 2003; Eisenhauer et al. 2003). The 
winds from these stars have a typical mass loss M w ~ 1O _6 M0 yr -1 and speed in the range 
v w = (1 — 3) x 10 3 kms^ 1 (Puls et al. 1996; Repolust, Puls & Herrero 2004). Loeb (2004) 
showed that SO-2, SO-16, SO-19 can supply the required mass flow that fuels the emission 
from Sgr A*. 

Now we can check if this mass accretion rate M w ~ 1O _6 M 0 yr -1 ~ 10 20 g s" 1 can 
provide the Eddington mass accretion into a SMMBH of ~ 50 g s -1 . In order to do this 
one should calculate the SMMBH luminosity of m smm bh = 3.7 x 10 -16 for a given color 
temperature T co i — 10 4 keV as follows (BRT99) 

Lsmmbh = 1.4 x 10 34 (7T 5 / 1 5) m 2 r [T co i (keV )/T h f ergs s -1 k2x 10 21 ergs s -1 (24) 

where r e // ~ 15 is the disk effective radius in Schwarzschild units, and the color factor 
T h ~ 2.6 (see BRT99). It is worth noting that we choose the color temperature T co i — 10 4 
keV based on the COMPTEL-EGRET observations (see the Introduction). 

Then we can evaluate the number of SMMBHs N sm mbh in the GC region if we assume 
they are responsible for 7— ray emission about 4 x 10 38 ergs s _1 , namely N smrn bh = 4 x 



12 - 


10 38 /2 x 10 21 = 2 x 10 17 . Note that the total mass of SMMBHs m to t = M to t/M Q inside the 
region of radius R ~ 3 x 10 14 cm is rn tot = m smm bhN smm t,h ~ 70. One can conclude that the 
mass flow rate M smm bh N smm bh = 10 19 g s~ 4 needed to fuel SMMBH in the Galactic Center 
region within 10 14 cm (or within ~ 10 3 Schwarzschild radii from the central BH) is much 
smaller then the mass flow 10 2 ° g s -1 supplied by the young massive O-stars in this central 
region. Consequently, we have shown how the small mass black holes could be become exposed 
when there is enough material around them. L 


6. Discussion 

6.1. The spectrum and origin of X-ray emission of a BH in the Galactic 

Center 

Revnivtsev et al. (2004) claimed that the inferred X-ray spectrum and luminosity of the 
emission illuminating the molecular cloud Sgr B2’s, strongly support the idea that the X-ray 
emission of Sgr B2 is Compton scattered and reprocessed radiation emitted in the past by 
the Sgr A* source. We have calculated the pair creation using the inferred luminosity and 
the inferred spectrum that is a power-law photon spectrum of index 1.8 with an exponential 
cutoff at 100 keV. In fact, this is a typical spectrum of BH sources at Iow/hard state which 
is presumably a result of Comptonization of the soft (disk) radiation in the hot electrons of 
temperature 50 keV in Compton cloud around central BH. 

It has already been shown that the spectral shape in the low hard state is independent of 
the luminosity (e.g. Laurent & Titarchuk 1999, Titarchuk & Fiorito 2004). We can propose 
that Sgr A* is now in a quiescent state, presumably a result of the gravitational energy 
release in the advection dominated flow (see e.g Narayan & Yi 1994). 

The luminosity of the X-ray spectrum of the X-ray source IGR J17475-2822 has not 
changed since GRANAT/ART-P and ASCA era, i.e. during last 10 years. The clear match 
of ASCA/GIS, GRANAT/ART-P and INTEGRAL/IBIS (Fig. 2, Revnivtsev et al. 2004) 
strongly supports this claim. Thus the proposed incident spectrum emanating from Sgr A* 
in the past (300 years ago) should not be vary either during the period at least of 10 years. 
It means that the inferred e+e- rate should be quite stable during the same period. The 
constraints on the “breathtaking variability on a time scale of only several hundred years” 
was proposed by Sunyaev et al. (1993). They argue that X-ray echo from Sgr B2 should 
be delayed by 300-400 years relative to the direct signal from Sgr A* due to the light travel 
time from Sgr A* to Sgr B2. One can put the important question regarding such a short 
recurrence time for the spectral transition, as 300-400 years for a black hole of mass ~ 3 x 10 6 



- 13 - 


solar masses. 

On the other hand one can suggest the more likely scenario for irradiating Sgr B2 was 
the interaction between the supernova remnant, we now see as Sgr A East and the 50 km 
s -1 molecular cloud behind Sgr A*, rather than Sgr A* itself. The luminosity, time scale, 
and spectrum all fit the requirements rather well, so the need for a high X-ray activity in 
Sgr A* can be weaker now. This effect of supenova (SN) ejecta and its interactions with 
molecular clouds proposed by Bykov (2002, 2003) is recently discussed in the literature (see 
e.g. Muno et al. 2004b, and Park et al. 2004). One can be right that “the luminosity, time 
scale, and spectrum all fit the requirements rather well” . However Park et al. pointed out 
that although the SN ejecta contribution to the observed neutral Fe line emission appears 
to be plausible, they noted a caveat, that nonthermal radio emission features have not been 
reported in the northeast regions of Sgr A*. 

This absence of nonthermal radio features is admittedly a difficult problem in terms of 
the SNR interpretation. It means that the supernova ejecta model needs the nonthermal 
distribution of electrons in order to explain the hard X-ray energy spectrum of the Sgr 
B2 molecular cloud (Revnivtsev et al. 2004). We want to emphasize once again that the 
spectrum seen from Sgr B2 is a typical spectrum of the low/hard state seen in many BH 
sources. It has already been proven long ago (see e.g. Sunyaev & Titarchuk 1980) that 
this hard spectrum is a result of thermal Comptonization of soft (disk) photons in the hot 
Compton cloud surrounding the central part of a BH. Thus one does not need any fine tuning 
or nonthermal electron distribution to produce such a spectrum, because in the outskirts of 
the GC BH there are plenty of molecular clouds, and we should see the effect of the reflection 
of the BH X-ray spectrum from these configurations. Thus the spectra of molecular clouds 
are naturally formed as a result of reflection of BH emission from molecular clouds. 


6.2. The site and the origin of the annihilation line 

In connection with the previous points, one can raise a question related to expected 
rapid variability in the e + e— annihilation line flux for such a compact BH region. However 
the observations over several decades seem to not show any variability at all. 

It is correct to say that the variability in the e + e— annihilation line flux would be 
quite high if this line was formed in the compact central region. However, we propose that 
the pairs are produced in this compact region and only then do they propagate through the 
Galaxy and the e+e- annihilation line is formed through the Galaxy. Thus variability of the 
e+e- annihilation line flux is related to the scale of a few kpc and it is not by chance that 
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observations over several decades seem to not show any variability at all. 

We infer that the ~MeV gamma ray region, where the pairs are created due to X-ray 
illuminations by the central source, is quite compact. Its radius is ~ 3 x 10 14 cm. This 
region (within 10 14 ~ 15 cm), is consistent with the COMPTEL map (Schonfelder et al. 2000, 
Strong et al. 1998). The full sky intensity map, longtitude and latitude profiles for 1-3, 3-10 
and 10-30 MeV show the strong peaks within a few degrees near GC. However, the spatial 
resolution of the COMPTEL map is not better than one degree and thus the compact region 
of ~ io 14-15 cm cannot be resolved in COMPTEL map. It is worth noting that Strong, 
Moskalenko & Reimer (2000) argue that there is an unresolved point-source population in 
the inner Galaxy and that this source is an important contribution to the emission around 
10 MeV. Furthermore, recently Strong et al. (2005) confirm that their model spectrum 
compatible with EGRET and Cosmic-ray measurements still lies about a factor 2 below 
the COMPTEL data, as found previously (Strong et al. 2000). There they propose that a 
contribution from compact sources is the most likely explanation. 

Producing ~ 10 42 s _1 pair flux (or more, depending on the duty cycle) has consequences 
beyond simply calculating the annihilation rate. The region surrounding Sgr A* is magne- 
tized, with magnetic field strength B of order 10 G inside 10 14 cm. Out to ~ 2 — 3 pc the 
field is turbulent (matching the turbulent plasma generated by wind- wind interactions from 
the Wolf-Rayet and O/B stars in this region). Melia & Falcke (2001) have suggested that 
the magnetic field strength of this region should be in the range of ^0.2 — 1 mG. 

In section 4.2 we review the positron propagation in the Galactic Center region and we 
find the positrons can pass through there. Consequently the sources of annihilation emission 
are presumably distributed over the Galaxy. 

Recently Liu, Petrosian & Melia (2004), hereafter LPM04, formulated and solved the 
propagation and acceleration problem in the turbulent magnetized plasma for a given pair 
flux Q , magnetic field strength B, plasma interaction time t p , density n, the size of the 
acceleration site R and particle distribution of the background plasma. They showed that 
acceleration of electrons by plasma wave turbulence in hot gas in the inner-most part of the 
black hole reasonably account for Sgr A*’s millimeter and shorter wavelength emission in 
the quiescent state and for the infrared and X-ray flares. In their model A (Q = 3.2 x 10 42 
s _1 , B ~ 9 G) they found the power carried away by accelerated electrons, with 7 > 100, is 
about 2 x 10 37 ergs s -1 ; which is more than enough to power observed radio emission, whose 
liminosity of 10 34 — 10 35 ergs s -1 . 

For the pair flux of order 10 42 s -1 , the “collateral” spectrum (and flux) arising from 
these particles is very consistent with Sgr A*’s millimeter and shorter wavelength emission 
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in the quiescent state and for the infrared and X-ray flares (LPM04). It is worth noting that 
the low energy electron flux Et me 2 calculated by LPM04, is comparable with the injected 
flux Q. It means that a significant fraction of the low enery electrons escapes from this 
region without any noticeable acceleration. In fact, the escape time of low energy electrons 
is more than a factor of 2 shorter than the time scales for acceleration and losses. The exact 
value of the electron flux depends on the model assumptions (see LPM04); and thus these 
theoretical uncertainties of the flux estimate (within a factor of a few) are unavoidable. The 
LPM04’s model also suggested that an outflow of high energy electrons (Liu & Melia 2002) 
are distributed over a spatially larger scale than that of the BH inner region. Finally the 
escaping high energy electron gas (cooled through the radiation) could form the annihilation 
line emission in the large scale of the Galaxy. The LPM04’s injected flux Q ~ 10 42 s -1 , is 
related to the mass accretion rate of order 10 -5 of the Eddington value (M Edd = L Edd /c 2 ). 
If the pairs are accelerated to 7 > 100 then their energy would be enough to power the 
observed radio emission whose luminosity is about 10 34 •*- 10 35 ergs s -1 . 

It is worth noting that within our model, the same flux of pairs (Q ~ 10 42 s _1 ) is 
generated in the SMMBH site as a result of X-ray illumination emitted by surrounding 
material. If at the present time the X-ray flux is of order 10 36 erg s _1 , then the pair flux 
would be 10 39 s -1 . 


6.3. The dark matter problem and its relation to observed annihilation line 

flux 

Dark matter is a long standing puzzle starting from Le Verrier with his discovery of 
Neptune in 1846 as a perturbation of Uranus (Le Verrier 1846). Our Galaxy is surrounded 
by an extended halo of unseen material. So far undetected, that dark matter induces a flat 
rotation curve in the Galactic plane (van Albada & Sancisi 1986). Its presence has been 
noticed since many decades (see Oort 1932). However, its nature is still unknown. The dark 
matter problem consists simply, in the existence of invisible mass showing its presence by 
gravitational effects. Now, it is widely accepted that dark matter exists. In fact, there is 
evidence for dark matter on scales from galaxy to the cluster of galaxies and to the whole 
Universe itself. Zwicky (1937) measured the dispersion velocity in the Coma cluster and 
found that the dynamical mass was hundred times more than the luminosity mass. 

Moreover, the dark matter seems essential: the growth of structure in the Universe by 
hierachical merging of the dark matter halo is a master piece of moderm cosmology. Recently, 
WMAP has established the presence of non baryonic dark matter with the density reported 
in terms of critical density of order Q m = p m /p c ~ 0.3. This is now the standard value of 
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the modern cosmology for dark matter density. 

On the galactic scale, one can interprete the rotation curve measurements by the pres- 
ence of the dark matter halo and then determine the mass density profile of this distribution. 
The average density is p m — 0.3 GeV cm -3 . The core of this halo seems to be extremely 
packed. In the Galactic Center, numerical simulations of structure formation in the non- 
linear regime have shown the presence of a singular power law cusp p(r) = /9o(r/7o) -7 , where 
0 < 7*2 and p 0 is the dark matter density at radius r 0 ~ 2 pc (Gnedin & Primack 2004), 
The central density p 0 is not uniquely determined and its value depends on models and it 
varies from 60 to 600 M 0 pc -3 . On the other hand up to now there is no identification for 
it despite the large number of candidates proposed from masses in the range between 10 -5 
eV to 10 12 M q . 

While there is no a priori reason to expect that massive compact halo objects (MACHOs) 
in the mass range 10~ 15 M Q <M < 10 -7 M 0 comprise a significant fraction of the density of 
the Universe, neither is there any definitive argument ruling them out. Nemiroff & Gould 
(1995) and Marani et al. (1999) discuss the possibility of identification of SMMBHs by 
gravitational lensing. The basic idea is to study the deflection of light by measuring the 
amplification due to gravitational bending of compact objects. For black hole mass in the 
range 10“ 16 M 0 to 10 -13 M 0 , the expected image separation is of order of 10 -15 arcseconds 
and the average time separation between the images is order of 10~ 24 s. They applied this 
delay as an interference process using gamma ray bursters. On the other hand, the absolute 
value of the light deflection by small black holes is extremely small (10 -12 cm)//?, where 
R is the characteristic distance in the Galaxy and therefore not observable. It is easy to 
show that the optical depth for the photons scattering at small angles by SMMBH is small, 
(~ 10 -6 ) and consequently the light multiple scattering by the small black holes is also ruled 
out. 


The supersymmetry theory as an extension of the Standard Model, predicts that one 
stable particle should be around today as a cosmological relic. Among these new particles, 
the neutralinos, could be the dark matter. If the dark matter halo of our Galaxy is made 
up of such particles, they will annihilate each other into fermion-antifermions pairs, and this 
will produce stable particles which could be detected as anti-protons, electrons and photons. 
Another way is to look for direct detection via laboratory experiments. Such new particles 
are scrutinized by direct and indirect searches since more than ten years. Our Galaxy is 
indeed a subject of intense activities because it could be a direct search for dark matter in 
studying astrophysical signatures. 

Among all the signatures, our Galaxy presents three strong signatures in X-ray, 7— ray 
backgounds and a typical e + e~ annihilation source in the Galactic Center. Among the 
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other possibilities, the micro-lensing searches have been developed and used to probe the 
dark matter halo of our galaxy (see above). The possibility that the halo could be made of 
primordial black holes is not ruled out by observations. Recently, Seto and Cooray (2004) 
discussed also the future possibility of searching dark matter candidates using the space- 
based gravitational wave detectors. They used this idea to test the possibility that the halo 
could be made of primordial black holes. According to Hawking predictions such primordial 
black holes will not suffer evaporation if their mass is greater than 10 15 g. 

While the consequences of evaporating black holes has been widely analyzed in terms 
of particles production, the case for stable SMMBHs s has not been considered so far. In 
literature, the problem is analyzed generically and the mass range used are often very large 
from 10 17 g to 10 20 g. We suggest the SMMBHs form the unseeing dark matter. 

7. Conclusions 

In this Paper, we show that there is a particular black hole mass which is very inter- 
esting for the dark matter problem. This mass of order of 10 17 g, corresponds to a typical 
Schwarzschild radius of order close to that of the proton. We explore this possibility in the 
light of existing galactic diffuse X and 7 backgrounds and with the strong source of annihi- 
lation lines at 511 keV in the Galactic Center. Our approach is to find a natural explanation 
of such astrophysical backgrounds within an astrophysics! dark matter scenario. 

This also can provide an explanation for origin of the annihilation line appearance in the 
our Galaxy which has been detected in several balloon and space experiments. We combine 
X-rays and 7-observations of the Galactic Center region to show that pair production can be 
a result of interactions of X-rays and gamma photons emitted in this region. Pair production 
strongly depends on the compactness of the X-rays and 7 emission region which is determined 
by X-rays and 7 luminosity and size of the region. 

For a velocity of the interstellar medium of much less than 10 8 cm s -1 and the value 
of density column a few times 10 22 cm -2 , the optical depth for 511 keV line is greater than 
one. We found for the observed X-rays and 7 luminosities and for the observed annihilation 
lines flux, the size of the X-ray and 7 emission area is ~ 3 x 10 14 cm. This size is consistent 
with the size of the area where the OB stars supply sufficient amount of fuel required for the 
X-rays and 7 emissions from the Galactic Center region. 

If we compute the total mass of SMMBHs inside the region of radius R = 10~ 4 pc, 
we obtain that this mass is about 70 solar masses or less (for this estimate see the end of 
section 5). One can then naturally ask if SMMBHs can be a candidate for dark matter 
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as already proposed in literature ( e.g. Nemiroff & Gould 1995, Marani et al. 1999). We 
compare this number to the dark matter mass computed using the singular power law cusp 
p(r) — Po( r A'o ) -7 (where 0 < 7<2 and p 0 — 60 — 600 is the dark matter density at 
radius r 0 ~ 2 pc) often used in the dark matter model for the Galactic Center (see Moore 
et al. 1998 and Gnedin & Primack 2004). It is evident that the value of the dark matter 
mass inferred from this model should be related to the radius of order a few parsecs and one 
should be careful to extrapolate this mass value to the radius of four orders of magnitudes 
less, namely to R == 10~ 4 pc. For this radius we find that the upper limit of the inferred 
mass using the dark matter model is at least of one order of magnitude less than that we 
obtain for SMMBHs. 

Koyama (1996) and Muno et al. (2004a), Revnivtsev et al. (2004) noted that the 
Galactic Center exhibited intermittent activities with time-averaged energy generation rate 
~ 1 q 41 ~ 42 er g S s -i comparable to Seyfert nuclei. Muno et al. (2004a) are left to conclude 
either that there is a significant shortcoming in understanding of mechanisms that heat the 
interstellar medium in Galactic Center region or that of a population of faint (< 10 31 ergs 
s _1 ) hard X-ray sources. Muno et al. (2004b) using a deep Chandra observations study 
to detail the origin of the point like sources near the Galactic Center. Possibly a single 
population of sources may dominate in this cluster. There are certainly many classes of 
objects present in this set of faint sources. We may suggest that some part of X-ray and 
7— radiation of the Galactic Center, can originate in accreting flows into SMMBHs (that are 
really faint sources with their bolo metric luminosity < 10 22 ergs s -1 ) in which outer parts of 
flow emit X-rays and the innermost part of the flow emit 7— rays. 

Finally we can conclude, that the only way to identify small mass black holes (SMMHHs) 
is by accretion mechanism in presence of a surrounding dense material. In the absence of this 
matter these SMMBHs are invisible and can be treated as dark matter, or invisible matter, 
i.e. without any observational appearance in the radiation. The crucial observational test of 
the SMMBHs as a model for a dark matter is detection of the blackbody radiation of ~ 10 
MeV temperature from the Galactic Center region. In the energy range between 1-30 MeV 
the photon energy flux of this emission should be of order 2 x 10~ 2 cm -2 sr _1 s _1 MeV which 
is compatible with that which has already been detected by COMPTEL. 
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Fig. 1. — Cartoon picture of compact emission area near the Galactic Center. Sources of 
' 7-ray region in blue, and sources of X-ray radiation in red. In upper and lower panels we 
show X and 7 rays spectra respectively. 





